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f^i ' Partial Split Supersymmetry with bilinear R-parity violation allows to reproduce all neutrino mass 

and mixing parameters. The viable dark matter candidate in this model is the gravitino. We study 
D , the hypothesis that both possibilities are true: Partial Split Supersymmetry explains neutrino physics 

and that dark matter is actually composed of gravitinos. Since the gravitino has a small but non-zero 
CNI ' decay probability, its decay products could be observed in astrophysical experiments. Combining 

bounds from astrophysical photon spectra with the bounds coming from the mass matrix in the neu- 
O I- trino sector we derive a stringent upper limit for the allowed gravitino mass. This mass limit is in 

Qh good agreement with the results of direct dark matter searches. 

PACS numbers: 14.60.Pq, 12.60. Jv, 95.35.+d 
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Split supersymmetry (SS) was originally proposed to address some of the most conspicuous problems of 
supersymmetric models, which are fast proton decay and excessive flavor changing neutral currents and CP 
violation [lfl. In SS the solution to these problems is accomplished by considering all squarks and sleptons 
very massive, with a mass scale m somewhere between the supersymmetric scale M susy and the Grand 
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Unification scale Mqut- One of the Higgs bosons remains light, as usual in supersymmetric models, as 
well as the gauginos and higgsinos, with all these particles having a mass accessible to the LHC M\. 

If R-Parity violation (RpV) is introduced in supersymmetric models, lepton number and/or baryon num- 
ber will be violated as well, inducing a potentially too fast proton decay [3]. Nevertheless, in SS the trilinear 
RpV couplings play little role in the phenomenology, with the exception being in the gluino decay rate. In 
this case, only bilinear RpV (BRpV) is relevant, opening up the possibility for a neutrino mass generation 
mechanism, without running into the danger of a too fast proton decay. As in any BRpV model, in SS-BRpV 
the atmospheric neutrino mass squared difference is generated by a low energy see-saw mechanism due to 
the mixing between neutrinos and neutralinos |4|,|5|]. Nevertheless, at one loop the only contributions to the 
neutrino mass matrix, coming from loops with the Higgs boson and neutralinos, is not enough to generate 
a solar neutrino mass squared difference |6, zl. Thus, an additional contribution to the model is needed |8] 
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or the model itself has to be generalized. 

In Partial Split Supers ymmetry (PSS) all squarks and sleptons have a mass of the order of the SS mass 
scale m, but both Higgs doublets remain with a mass at the electroweak scale 17|,|9J]. The addition of RpV to 
this model was introduced to be able to generate a solar neutrino mass 170 . Loop contributions from neutral 
CP-even and CP-odd Higgs bosons are indeed able to do the job, producing not only the atmospheric and 



loo. 



solar masses, but also the atmospheric, solar, and reactor neutrino mixing angles 

Since R-Parity is not conserved, the lightest supersymmetric particle (LSP) will be unstable and decay 
into SM particles. This will occur very fast if the LSP is the traditional neutralino, losing it as a candidate 
for Dark Matter. It is however known that in the case of a gravitino LSP, despite being unstable, it will 



12|]. Suppression of the gravitino interactions 



decay very slowly via gravitationally induced couplings 111 ll . 
by both the Planck mass and the R-parity violating couplings leads to a very long-lived massive particle, 
whose lifetime can typically be of several orders of magnitude longer than the age of the universe. The next- 
to-lightest superparticle (NLSP), on the other hand, has a lifetime which is much shorter than 1 second, and 



thus becomes harmless to a successful big-bang nucleosynthesis IU3I1 . 

The fact that the gravitino decays to ordinary particles in such a scenario has given rise to interesting 
phenomenological studies jl4l4l6ll. An important one is the study of the photon spectrum produced from the 



two-body decay QJ/7 



19], and more recently from the three-body decays of the gravitino [20 



2l|]. Important 



constraints on the mass and lifetime of the gravitino can be derived from the fact that its decay has not been 
detected by gamma-ray telescopes, most importantly the Fermi Large Area Telescope (Fermi LAT) 1 22 -I26 1 . 

Further constraints on the allowed parameters can be derived from the measurements of the neutrino 
masses and mixing parameters, which in PSS provide constraints on the R-parity violating couplings 

In this work we study these two independent constraints, and show that when taken together they imply 
the existence of a maximal value of the light gravitino mass m^f, an d also the existence of a minimum 



value of the gravitino lifetime. Our findings are in good agreement with direct dark matter searches Q27 



28] 



which put much stronger constraints on dark matter particles with masses 3> mylf- It is important to men- 
tion that this kind of model can be further studied in the context of the observed Pamela electron/positron 



excess and direct LHC signals Q2S 



Further studies are possible in the context of the early universe 



1 3114320 . However, since in our work the gravitino has an extremely long lifetime, the implications from 
gravitino decay in the early universe are not taken into account. 

This model is a bottom up approach, it can however be motivated by a number of general considerations. 
In ref. it is shown that a Split Supersymmetric spectrum can be easily generated in models with direct 
mediation of supersymmetry breaking, i.e., without invoking a hidden sector. In opposition to low energy 
supersymmetric models, in SS this is possible because the gauginos are allowed to be much lighter than the 
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sfermions. A toy model is given, where sfermion masses are generated at a scale fh while gaugino masses 
are induced at a lower scale //, after integrating out physics at a higher scale that controls the supersymmetry 
breaking. In the same model, a gravitino mass is generated at the scale fh 2 /Mpi allowing it to be even lighter 
than gauginos. This toy model generates also a term of the order of fh 2 and it would have to be modified 
in order to accommodate PSS, where is required to be much smaller than fh 2 . It is worth to mention 
that a B^ <C fh 2 would alleviate the otherwise present fine-tuning needed in the Higgs potential to generate 



a correct electroweak symmetiy braking 113511 . As a top down approach to PSS we mention also ref. |9|] 
where a PSS spectrum is given. In this model sfermions masses are generated at a scale fh 2 ~ Vq/M 2 , 
where supersymmetiy breaking originates with a hidden sector dynamics with vacuum energy Vq, and it is 
communicated to the visible sector by massive fields of mass M. In addition, two light Higgs scalars are 
composite with a mass ~ ■£-A. comp , with A cornp the typical scale of the strong dynamics. The gravitino 
mass is of the order \/Vq/Mpi and can be the LSP for values for example Vq ~ (10 9 GeV) 4 . 



II. PARTIAL SPLIT SUPERSYMMETRY 



In PSS both Higgs doublets remain with a mass at the electroweak scale. As it happens in SS, higgsinos, 
gauginos, and Higgs bosons interact via induced couplings of the type, 

JC% C S 3 - j=Hl(g u aW + ~g' u B)H u - j=H j d (g d aW - ~g' d B)H d + h.c, (1) 

where g u , g' u , g d , and g' d are couplings induced in the effective low energy lagrangian. At the SS scale fh 
they satisfy the boundary conditions, 

(2) 



9u =9d = 9- 



9u = 9d = 9, 



evolving with independent RGE down to the electroweak scale. Similarly to the MSSM, both Higgs fields 
acquire a vacuum expectation value (H u ) = v u and (Hd) = Vd, with the constraint v 2 = v 2 + v 2 , = 
246 GeV 2 and the definition tan/3 = v u /vd- Gauginos and higgsinos mix forming the neutralinos, with a 
mass matrix very similar to the one in the MSSM, 
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(3) 



where M\ and M 2 are the gaugino masses associated to the U(l) and SU(2) gauge bosons, and fi is the 
higgsino mass. In our calculations we will neglect the running of the higgsino-higgs-gaugino couplings and 
work with their approximated value indicated by eq. ©. 
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The addition of R-Parity violation to the PSS lagrangian allows us to study a mechanism for neutrino 
mass generation. Trilinear couplings are not relevant for this problem, because all squarks and sleptons are 
very heavy, with a mass of the order of fh, and thus decoupled from the low energy effective theory. Only 
BRpV is relevant, and is described by the terms, 

Zpfs = -dH^Li - ^biHle{g d aW - ~g' d B)U + h.c, (4) 

where e = iai- The are the supersymmetric BRpV parameters in the superpotential, which at the low 
scale manifest themselves as mixing between higgsinos and lepton fields. The 6j are three dimensionless 
parameters attached to lepton-higgs-gaugino interactions. They are analogous to the ones in eq. (Q]), except 
that violate R-Parity, and are generated in the effective low energy theory. 

The origin of the BRpV terms in dU is related to the /i-problem which refers to the origin of the term 



[iH u Hd in the superpotential. As discussed for example in ref . IB6I1 , the same mechanism that solve the 
//-problem could be used to explain the origin of the e terms. A popular mechanism is the existence of non 
renormalizable couplings of the sort (a/Mpi)&i&2H u H c i, with a being a dimensional constant and <E>i and 
$2 two hidden sector scalars. When supersymmetry breaks and <&i and $2 acquire vacuum expectation val- 
ues a //-term is generated. Similar terms (ai/Mpi)^^iH u Li can be present generating e-terms, although 
di and/or ($3), ($4) should be much smaller than a and/or ($1), ($2) in order to have €, L a few order of 
magnitude smaller than fi, necessary for neutrino physics. The b{ terms are generated after integrating out 
the sleptons, and appear because above the scale m the Higgs bosons mix with sleptons. At the scale m we 
have bi ~ Vi/v u , with the necessary condition v u B e i ~ ViM^ from the minimization of the scalar potential 
fl. Large B ei can be easily obtained in a similar way as for B„, as explained in fl. We mention also 



ref. 0711 where R-Parity is naturaly broken radiatively when a right-handed sneutrino acquires a vacuum 
expectation value, generating bilinear R-Parity violating terms. In our work, we concentrate on the effect 
of bilinear R-Parity violation. Trilinear R-Parity violating couplings could be present, but we ignore their 
effects due to the large mass of the sfermions. 

III. NEUTRINO MASSES 

When Higgs bosons acquire vacuum expectation values, mixing terms between gauginos and leptons 
are generated, producing the following mixing between neutrinos and heavier fermions, 



r RpV 
'-PSS 



EiHl + -hv u (g d Wz - g' d B 



Vi + h.c. + . . . (5) 
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The effect of these terms is that neutralinos mix with neutrinos, forming a 7 x 7 mass matrix, which in the 
basis ( B, W, H u , Hd, v e , v T J has the form, 



M 



N 



(6) 



m 

and where the 4x4 submatrix M x o is the neutralino mass matrix in eq. ©. The 4x3 neutralino-neutrino 
mixing block 



m 



\g d biv u \gdhv u ei 



(7) 



develops from terms in eq. (|5). A diagonalization by blocks of the mass matrix in eq. © can be achieved 
by a rotation AfM. at AA T given by, 

-Z 1 



(8) 



where we define, 



with 



£ = mM- \ 



(9) 



&4 



2 det M x o ' 



A; 



3d pM: 



l -A< 



2 det M x o 



(10) 



■ (Miq u q d + M 2 q'q'd) K - — 
4detM x o v iyuyd ^ ^ 



Vd 



4 det M x o 



{M 1 ~g 2 d + M 2 ~g'i)K. 



and Aj = \ibiV u + e^. This leaves an induced effective neutrino mass matrix equal to, 



m(°)|.- 



with 



A (0) 



4 det M x o 



(ID 



(12) 



At this level only one neutrino acquires mass, leaving the solar squared mass difference null and the solar 
angle undetermined. Quantum collections contribute to the neutrino mass matrix in such a way that the 
degeneracy in eq. (fTTT> is lifted, leaving it with the following form, 



M 



AAiAj +Ce i e i , 



(13) 
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where the B term [10] has been made to vanish by an appropriate choice of the subtraction point. Thus 
at one loop two neutrinos acquire a mass while the third one remains massless. Since in this case the 
experimental value Am^/Am^, ~ 0.035 implies m v . A » m U2 we have, 

2 



Am 



atva 



Am: 



ol 



A\A\ Z + C\e 
A 2 C 2 \A x e | 4 
(a\K\ 2 + C\Z\' 



2AC\K x e\ z , 



(14) 



For later use, we introduce the photino 7 and zino Z fields by rotating by the weak mixing angle the 
weakly interacting gauginos B, W, in direct analogy to their standard model counter parts 



/ 7 \ 

Z° 



I 



cw sw ■ ■ ■ 

-sw cw ■ ■ ■ 



/5> 
w 



J 



A 



w 



( b\ 

w 



(15) 



where the dots indicate that all other states are just multiplied by the unit matrix. Thus when dealing with 
this new basis (7, Z°, H u , H^, v e , v^, v T ) the mixing matrix is 



M'=MAl 



w ) 



(16) 



where only the first two states are rotated. 



IV. GRAVITINO DECAY 



Indirect observation of the gravitino becomes a possibility due to its decay to ordinary particles. In this 
section we calculate the possible decay channels of the gravitino as dark matter candidate, assuming that 
m 3 / 2 < mw- We then relate these results to experimental bounds on the decay products and on the RpV 
parameters. 



A. Two-body decay 

When R-parity is conserved the gravitino can radiatively decay into a photon via the following term in 
the Lagrangian, 

where Mp is the Planck mass, ifi^ is the spin-3/2 gravitino field, A 7 is the spin-1/2 photino field, F vp = 
d v A p — dpA u is the photon field strength, and A^ is the photon field. This coupling might in principle 
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be modified by a factor of order one. In this work however we assume that (fTTT ) is the exact form of the 
coupling. Variations of the final result can then be studied by order one shifts of the Planck mass Mp. This 
Lagrangian term induces the following decay, 




where the cross indicates we are picking the photino component of the corresponding neutrino, which mix 
due to violation of R-Parity. The amplitude for this decay is 

where U~ v is the amount of photino in the neutrino fields, as indicated by the neutrino eigenvector. We 
write the gravitino field as the tensor product of a spin- 1/2 field with a spin-1 field, 

%(p, A) = £(1/2, s; 1, m|3/2, X)u(p, s)e lt (p, m), (19) 

s,m 

obtaining the following completeness relation jn|, 



_ PfiPv \ _ 1 / _ PfiPa \ ( _ PvP\ \ a A 

9flU m l/ 2 J 3 \ m 3/ 2 /V m l /2 P 1 , 



■ (20) 



A 

Knowing the above relation, the calculation of the differential cross section is standard, giving the result 

(\M \ 2 ) ( m 3 /2\ _ m\ /2 / m 3/2 



dE^dSl ~ 64^m 3/2 6 \ Ej 2 ) ~ 128tt 2 MJ, ^ H^ 7 2 J' (21) 
which is independent of the angles as expected. The total decay rate is then 

T(G^ju) = ^V|^J 2 . (22) 

For a gravitino mass m 3 / 2 < mw this is the only kinematically allowed 2-body decay. By using the relation 
in eq. ( fT6l ) one finds that the photino-neutrino mixing factor in PSS is 

U^ v . = Nucw + Ni2S\v, (23) 
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with i = 5, 6, 7 labeling the neutrino generation, and where tw = (f /g is the tangent of the weak mixing 
angle. Using equations © to (fTOl) we find, 

U ™ ~ 2 (detM o) ^ dMlSw ~ 9dM2C W ) K- (24) 

For the numerical calculations we will require the sum over the generations of the square of the mixing 
factor: 

2 

\U^f := J2 4* " 4 ( det M x0)2 (9 d M lSw - ~g' d M 2 c w f |A| 2 . (25) 
This is because we do not distinguish the different neutrino flavors. 



B. Three-body decays 

When studying the three body decay a more general part of the interaction Lagrangian comes into play 

C 3 [(D^n^YPLX* ~ [D^)x 1 Pr1 v 1^] (26) 



where the second line is in analogy to (flTT ) and the first line introduces additional couplings with scalar 



2111 . where 



fields cp. The 3-body decays of the gravitino were studied in detail for the first time in (20 
explicit formulae are given. Nevertheless, our calculations have yielded that the three-body results in ||20H 
have to be corrected. We agree, however, with the conclusion that the 3-body decays are indeed important, 
and cannot be neglected. We find 3-body decay branching ratios of the order of 10% for gravitino masses 
of order 10 GeV, and greater for larger masses. The exact formulas for the amplitudes of the contributing 
diagrams are given in the Appendix. 

First we consider the gravitino decay into a fermion pair and a neutrino. The first pair of diagrams are 
three-body decays via an intermediate photon and Z boson, 
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where the cross means we take the photino (or zino) component in the neutrino field. These two amplitudes 
are equal to, 



Ml = i^k 1 )(-iq f h,v(k 2 M q )^^r(h' 1 -7 ,i m a (p), 



(27) 



M 2 = i 



(k z — m z ) + imzTz cyy ^ iMp 



with Cy = Tg /2 — qfSyy and c A = —Tl/2. Their contribution to the decay rate is given in the Appendix. 
The photino-neutrino mixing is given in eq. (I24l ). while an analogous calculation for the zino-neutrino 
mixing gives, 



(g d Mic w + g' d M 2 s w ) Ai 



(28) 



Zv - 2(detM x o) 

To the previous two amplitudes we add a contribution coming from quartic couplings between the gravitino, 
a gauge boson, and a scalar with its fermionic partner, with the scalar acquiring a vacuum expectation value, 




/ + 




i Mi 



The first amplitude is proportional to b{V u and the second to v^fe/ fi), in such a way that the combined 
amplitude Ais is proportional to Aj, 



M 3 = 



[k z — m z z ) + imzl z cw AcwfJ-Mp 
whose contribution to the decay rate is also given in the Appendix. Since the neutrino is not directly 
detected, a sum over flavors must be done [as in eq. (l25lll. after which it is clear that the decay rate will 
satisfy T(G -> ffv) oc |A| 2 . 

Now we consider the gravitino decay into two fermions and a charged lepton. The decay via a W gauge 
boson is represented by the Feynman diagram 
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whose amplitude can be shown to be, 

M A = 



—i 



%(1 - 75 ) v (A ;2 )n( (? )^5^7 a (^ " l^mUp)- 



(k 2 - m 2 w ) + im w Y w v iy 2y/2"*'~ ""'v*' ^ 4M P 
Its contribution to the decay rate is given in the Appendix, and it is proportional to the wino mixing to 
charged leptons, given by, 

9d 



-A, 



(30) 



wti V2{detM x+J 

where M x + is the chargino mass matrix [7]. This graph is complemented by an amplitude coming from 
quartic couplings between the gravitino, a W gauge boson, and a neutral scalar with its charged fermionic 
partner, with the scalar acquiring a vacuum expectation value, 





iMt 



As before, the first amplitude is proportional to b,iv u and the second to Vd(ei/fi), such that the combined 
amplitude A^ 5 is proportional to Aj, 

—i 



Mr = i 



(fc 2 - m 2 w ) +im w T w "' v " / 2V2 '°'" v ""^' 2^M P 



11 



1 

0.8 
* 0.6 

O 

I 

S 0.4 
0.2 



" 20 40 60 80 

m 3/2 [GeV] 

FIG. 1: Branching ratio of the three-body decay of the gravitino as a function of its mass TO3/2, red (full 
line) for M\ = 100 GeV, green (dotted line) for M\ = 300 GeV, and blue (dashed line) for M\ = 500 GeV. 

Its contribution to the decay rate is also given in the Appendix. Note that when fermion / is a charged 
lepton and /' is a neutrino, we get interference between diagrams .Mi 2,3 i n one hand and in the 

other. This is because the decay G — > i'ftTVj can proceed via a Z or a W gauge boson. In addition, note 
that r(G -»• //\~) = Y{G -> /'/^ + ), thus we multiply the first one by two. 

The decay rate for the 3-body decays of a gravitino can be written as a sum of the various terms given 
in the Appendix. We sum over three generations, and neglect the masses of the final states. The end result 
is proportional to |A| 2 . Thus, in the branching ratio the factor |A| 2 cancels, depending only on M\, M2, 
and the gravitino mass m 3 / 2 ifll - Here we work with the assumption Mi ~ 2M\. The result is shown in 
Fig. [J where we plot the branching ratio as a function of m 3 / 2 for the three values of Mi = 100, 300, and 
500 GeV. One observes that the three-body decay becomes important for large gravitino masses m 3 / 2 and 
large M±. The dependency on m 3 / 2 can be understood as a phase-space effect in the three-body decay rate. 
The influence of Mi on the branching ratio can be understood by the fact that the two-body and three-body 
decays, by the virtue of the mixings Uj V and U% v respectively, get suppressed by growing M±. In contrast 
the parts of the three-body decays that contain vacuum expectation values (z^) and (Hd) do not experience 
this suppression and thus become more important in the regime of large M\. Those effects are also present 
in it is only the form of the curves that turned out to change in the corrected version. For the masses we 
are interested here, 7713/2 ~ 10 GeV, the 3-body decay is ^$ 10% and with small dependence on the gaugino 
mass. For larger gravitino masses the branching ratio can be as large as 80%. Since the calculations are 
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being carried out in the Feynman gauge, there are in principle also diagrams containing Goldstone bosons 
in the propagator. However, due to the coupling of the Goldstone bosons, those contributions vanish in the 
limit of light fermion masses. 



C. Induced photon flux 



The photon spectrum produced by the decay of the gravitino consists of a mono-energetic line of energy 
777,3/2/2 from the two-body decay, plus a continuum distribution from the three-body decays. The exact 



form of the spectrum, which depends on m 3 / 2 and Mi, was studied in detail in |20j, 



2111 using an event 



generator. Here we are interested in obtaining constraints on the gravitino parameters, for which it suffices 
as an approximation to consider only the photon line from the two-body decay, as this is the most prominent 



feature of the spectrum for values of Mi up to 1 TeV [20]. Including the three-body decay would make our 
final conclusions slightly more stringent. 

The observed spectrum is calculated from the flux of gamma-rays expected at earth. This flux is the 
sum of two contributions, one from the gravitinos decaying in the galactic halo, and one from the gravitinos 
decaying at cosmological distances. It has been shown that the first contribution is highly dominant, and so 
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we will neglect the second one 



has the following simple form [[18 



. In this way the differential flux, as a function of the photon energy E, 



19]: 



E 



2 dJhaXo 



D, 



7 m 3/2 



dE 



5 ( E 



m 3 / 2 



(32) 



with 



D- 



r(G -> 71/) 



Phaio(r)d£ ) = d 1 T{G ->■ jv). 



(33) 



8vr Vl.o.s. 

The constant d 1 depends on the dark matter density profile of the halo, and on the region of the sky consid- 
ered for averaging the flux (denoted by the term in brackets above, where l.o.s. means line of sight). Using 
the Navarro-Frenk- White profile 113 811 . and considering the region \b\ > 10° for the average (with b denoting 
the latitude in galactic coordinates), we find 



d 7 = 0.80 x 10 24 [MeV cm -2 str -1 ]. 



(34) 



This region of the sky was the one considered by the Fermi LAT collaboration in the derivation of the 



extragalactic diffuse spectrum [23]. 
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V. LIMITS AND CONSTRAINTS 



A. Constraint from the observed photon spectrum 



The fact that the line produced by the gravitino two-body decay has not been observed gives constraints 
on the mass and the lifetime of the gravitino. Assuming that the extragalactic diffuse spectrum measured 



by Fermi LAT can be correctly modeled in terms of known sources j22ll . one can use this spectrum to find 



constraints 0231 



After the convolution between the calculated flux in eq. (1321 and a Gaussian distribution we find 

E 2 dJ telo = D l m 3/2 1 -{E-m 3/2 /2?/2^ (35) 
dE 2 

where a is related to the energy-dependent resolution of the Fermi LAT instrument p evaluated at the 2-body 
peak 

m M2 

a = pE = p^- . (36) 



The energy dependence of p evaluated at E = m 3 / 2 /2 can be approximated by 124 



3i 



p(m 3/2 ) = 0.349 - 0.142 log + 0.019 log 2 (^L) . (37) 



Thus, the maximum of the photon spectrum is given by 



-,2 <^halo \ 



# — 7^ = — j= (38) 

/ max 

P( m 3/2)V2vr 

On the other hand, the intensity or the integrated flux for the extragalactic diffuse emission (fo r the range 



E > 100 MeV and the sky region |6| > 10°) was measured by the Fermi LAT collaboration B23I1 . 

f°° dJ 

/(> 100 MeV) = / — = (1.03 ± 0.17) x 10" 5 [cm" 2 s" 1 str" 1 ]. (39) 
Jwo dE 

together with the observation that the spectrum can be fitted by a power law dJ/dE oc E~~* ', with index 
7 = 2.41 ± 0.05. From this we calculate the spectrum to be 

p2 dJ tac I i c - ! „ 1 3 fMr^r „ m — 2 —1 ll ( ^ \ 



E z — = (9.6 ± 1.6) x 10~ d [MeV cm" 2 s" 1 str" 1 ] — — , (40) 
dE y ' 1 J I 1 MeV J 



where the error was directly calculated from the error in the integrated flux using eq. (1391 . 

As we mentioned before, we assume that the central value of the spectrum can be explained by models 



of known sources [22]. We impose that the extra contribution from the gravitino source is smaller than a 3a 



error margin. This gravitino contribution is related to the decay rate T(G — > ju) through eqs. (|33l and (|38l ). 
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By introducing the total decay width, r 3 h = T(G — > 71/) + T(3 — body), we find the following restriction 
on the gravitino lifetime, 

r3/M >4^(G^)^r 2 , (4i) 



10 27 s/ p(m 3/2 ) v ' ; VI GeV 

where B(G — > 71/) denotes the branching ratio of the 24)ody decay. In this way eq. (|4TI) defines a region in 
the m 3 / 2 — T3/2 plane consistent with the non-observation of the gravitino decay by Fermi LAT. 

The gravitino 24jody decay width F(G — > jv) is given in eq. (l22l) . If we sum over all neutrino species 
it becomes proportional to |A| 2 , via |?7^| 2 in eq. (T25T ). The 34x>dy decay width is also proportional to |A| 2 . 
The reasons are analogous to the 24jody decay, since after summing over lepton generations we see that 
|.Mi| 2 is proportional to |£7^t,| 2 , |.M 2 | 2 is proportional to \U^J 2 [eq. (1281)1. and |A^4| 2 is proportional to 
\U^ £ \ 2 [eq. (l30l)l. In addition, amplitudes |.M 3 | 2 and | 2 are directly proportional to |A| 2 as can be 
seen from eqs. (1291 ) and (|3TI ). In this way, the gravitino lifetime becomes large for two reasons, because the 
Planck mass is large and because BRpV is small: t~^ 2 oc |A| 2 /Mj,. 

We display experimental constraints on the model in the m 3 / 2 — t 3 / 2 plane, with the first one given 
by eq. (|4TT) . This constraint depends also on \U^,y\ 2 and the gaugino masses. In fact, the whole decay 
rate T(3 — body) can be factored out by \U~ /Iy \ 2 with the remaining factors depending on M\ and M 2 , but 
with the dependence on \i being in first approximation negligible. We further use the simplifying assumption 
M 2 = 2M\ motivated by mSUGRA models. In this way, we display the constraints in the plane m 3 / 2 — r 3 / 2 
as a function of | U~ v | and Mi . In the first constraint in eq. (|4T1) though, the dependence on | U~ v \ drops out. 

B. Constraints from the Neutrino Mass Matrix 

Further constraints appear from neutrino physics, controlled by the BRpV parameters Aj and ej, and by 
MSSM parameters like gaugino and higgsino masses. We do a scan over parameter space looking for good 
solutions to neutrino observables. The range in which we vary the parameters is given in Table U 



We define a % 2 value for each point in parameter space as follows 



X 



/ 10 3 Am 2 tm -2.4 
V 0.4 



+ 



10 5 Am 2 ol - 7.7 
06 



+ 



sm 



0.505 \ 2 



0.165 



sin 2 e soi - 0.33 

0T07 



(42) 



allowing a 3<r deviation 14 111 . The point is accepted if \ 2 < 4, plus the additional condition that the reactor 



angle satisfies the bound sin 2 # reac < 0.05. Since r 3 / 2 depends directly on \ U^ U \ we deteimined its maximal 
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SUSY parameter 


Scanned range 


Units 


tan/3 


[2,50] 


— 


iMl 


[0, 1000] 


GeV 


M 2 


2Ml 


GeV 


Mi 


100,300,500 


GeV 


m h 


[114,140] 


GeV 


m A 


[50,6000] 


GeV 


Q 


951.7 




RpV parameter 


ei 


[-1,1] 


GeV 


£2 


[-1,1] 


GeV 


£3 


[-1,1] 


GeV 


Ai 


[-1,1] 


GeV 2 


A 2 


[-1,1] 


GeV 2 


A 3 


[-1,1] 


GeV 2 



TABLE I: Scanned ranges for PSS and RpV parameters. 



and minimal values for a given Mi compatible with neutrino physics. The numerical results are given in 
Table HI 



Mi 


|[/^| 2 (min) 


|?7^| 2 (max) 


100 GeV 


2 x 10~ 16 


4 x 10~ 13 


300 GeV 


2 x 10~ 17 


3 x 10~ 14 


500 GeV 


1 x 10- 17 


1 x 10- 14 



TABLE II: Maximal and minimal values of |C/^^| 2 , consistent with neutrino experiments, for three 

different values of Mi. 



One sees that the range of possible values for \U^ U \ 2 depends on the value of Mi, with the maximal value 
being around 3 orders of magnitude greater than the minimum value for each case. Since the gravitino 
lifetime t~^ 2 = T(G — > ju) + T(3 — body) depends on [JT^I 2 , m 3 / 2 , and Mi, this imposes two extra 
constraints in the plane m 3/ / 2 — r 3 / 2 that complement the one in eq. (|41~T) . 
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VI. COMBINED CONSTRAINTS AND RESULTS 

The combination of the constraints found in the previous section defines an allowed region in the m 3 / 2 — 
r 3 / 2 plane. This region is shown in figs. (f2ab and (l2bl for the gaugino mass, Mi = 100, and 500 GeV. One 
sees that the constraint from the photon spectrum, when taken together with the maximal value of | U^ v \ 2 
consistent with neutrino experiments, gives a lower bound on the gravitino lifetime. In all cases that we 
studied, this bound is several orders of magnitude larger than the age of the universe, compatible with a 
good candidate for dark matter. Even more interestingly, we see from these graphs that the constraint from 
the photon spectrum, when combined with the minimum allowed value of |^i/| 2 , imposes an upper bound 
on the gravitino mass. This bound is near 2, 4, and 5 GeV for Mi = 100, 300, and 500 GeV. 

We also note from these graphs how the 3-body decays of the gravitino become more important as M\ 
increases. In particular, the constraint coming from the photon spectrum analysis becomes less stringent 
as the gravitino mass gets closer to myv, which is quite evident for the case of M\ = 500 GeV. This is 
expected from the fact that the strength of the gravitino photon line is proportional to the branching ratio 
of the 2-body decay. Conversely, for a gravitino of mass below ~10 GeV this branching ratio is close to 1, 
and so from eq. (l22l we see that the lifetime is approximately proportional to m%j 2 , as can be noted in the 
graphs. 




m 3/2 [GeV] m 3/a [GeV] 



(a) Allowed region for Mi = 100 GeV. (b) Allowed region for Mi = 500 GeV. 

FIG. 2: Allowed (shaded) region in the m, 3 / 2 — T3/2 plane. The region above the nearly horizontal (red) 
line is allowed by the constraint in eq. rt47D . The region between the oblique lines is allowed by minimum 

(green) and maximum (blue) values of\U^ v \ 2 . 

We stress that this analysis assumes that m 3 / 2 < myv- For a gravitino with a mass greater than the 
Z boson mass the two-body decays G — > Zv and G — > Wi will be kinematically allowed. As shown in 
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ref. Ill 51 . the branching ratio of the decay G — > becomes very small for a gravitino mass above 100 GeV, 
and so the monochromatic line in the photon spectrum becomes less important. A more detailed analysis 
of the photon spectrum produced by the decay of the gravitino would be required in that case, which would 
need to include the fragmentation of the W and Z bosons in addition to the contribution from the 3-body 
decays. 



VII. SUMMARY 



It is explained how PSS in combination with RpV allows to generate the neutrino masses and mixings 
at the one loop level. Then it is investigated whether in this model the gravitino still is a good dark matter 
candidate. In order to do this the gravitino decay rates into two and three-body states are calculated. In 
the three-body decay corrections are found to previous calculations. Since for relatively small gravitino 
masses the decay is dominated by the two-body decay, the corrections in the three-body case do not affect 
the final result. Those rates allow to calculate the additional photon flux that is induced by the gravitino 
decay. Comparison of this photon flux with recent data from the Fermi LAT collaboration, allows to put 
restrictions on our model in the gravitino decay process. 

Finally, combining the restrictions obtained in the neutrino sector with the restrictions in the gravitino 
dark matter sector of the same model, an upper limit on the gravitino mass rra 3 / 2 < 6 GeV is found. The 
exact values of the maximal m 3 / 2 and the minimal r 3 / 2 in our model are given in table JII] One observes a 
relatively weak Mi dependence. 



M 1 [GeV] 


r 3/2 (min) [s] 


m 3 / 2 (max) [GeV] 


100 


4.7 x 10 28 


2.3 


300 


7.2 x 10 28 


4.4 


500 


8.6 x 10 28 


5.3 



TABLE III: Minimal values of r 3 / 2 and maximal values of m 3 / 2 as a function of M\. 



It is interesting to note that most direct dark matter search experiments disfavor the typical heavy dark 



28fl. The bound, however, 



matter particle that appears in R-Parity conserving supers ymmetric models (27, 
is not even closely as tight for a light dark matter gravitino as it is found here. On the other hand, with 
the derived limit, the model turns out to be directly testable. For instance, if a relatively heavy dark matter 



particle is found with m 3 / 2 3> 10 GeV, one could immediately conclude that this model is ruled out. This 
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feature of testability and the possibility of falsification can be seen as a strong advantage over many other 
models. 
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VIII. APPENDIX: THREE-BODY DECAYS FORMULAS 

We label the five relevant diagrams with indices 1 to 5 as indicated in the text. For the photon and Z 
mediated diagrams we define the invariant masses s := (k\ + /C2) 2 , and t := {k\ + q) 2 , where k±, k 2 and 
q are the 4-momenta of the fermion, antifermion and neutrino, respectively. For the spin-averaged squared 
amplitudes and interferences we find (we define m = m 3 / 2 for simplicity), 



1 ( q)\U^\ 2 \ l_ 
16Ml ) s 2 



^=\(^k) (s - m |)Wi T22 ' (44) 



(l-^sl 2 ) 4 (^ 64/x 2 c / M | ) (s - m|) 2 + m|r| T33 ' (45) 

2K{ Ml M 2 ) - - y i6cwM , ) s[(s _ m|)2 + m|r|] ^ (46) 

osp/ U * M \ - 2 ( g 2 g/^ A A (f ~ ™ 2 Z ) 

2^ Ml Ms) - - 4 {-^r^r) s[(s _ m | )2 + m|r|] r u , (47) 

2 / g 3 U? Ai \ i 

2K(M^ 3 ) = 7 7^777 7 o75 ?7^ T 23- (48) 

For the W-mediated diagrams, we define s and t as above, with fci, k 2 and g the 4-momenta of the 
neutrino, antilepton and lepton, respectively. We find, 



19 



(| - M5|2) " \ (25^) (s-m*J + m^ w ^ (50) 

Finally, for the interference terms between these two groups of diagrams (when f = £' = £), we define 
s and t as above, with k\, k 2 and q the 4-momenta of the lepton, antilepton and neutrino, respectively. We 
find, 

2«<A< 1 A<4> - z ^ 32V - 2M2p ) a[(m 2_ (a + t) _ m ^ )2 + m 2^ ] ^ (52) 

?si?/A/p yu \ = - 1 92u zvi U wti \ (g - "^l)(^ 2 - (g + *) - ^vf) + m z miyr z r w 

^ 2 U ' 4 [32V2c w M 2 P ) [(s - m|) 2 + m%T%) [{m 2 -(s + t)- m 2 w ) 2 + m 2 w T 2 w ] 24 ' 

(53) 

?mM*M \ = - I g3 ^j Ai 1 (g - m|)(m 2 - (a + t) - m 2 ^) + TO Z mn/rzr w 

V 3 4/ 4 ^G^c^M 2 J [( S - m|) 2 + m|r|] [(m 2 - (s + t) - m 2 w ) 2 + m 2 w Y 2 w ] 34 ' 

(54) 

2R {Ml M 5 ) - - [ r ^ r ) s[{ ^- {s + t)-m 2 w ) 2 + m 2 w T 2 w ] Tl5 ' (55) 

WtM*M \ = -l '^EzuAA {s-m 2 z ){m 2 -{s + t)-m 2 w )+m z m w T z Y w 

K{J hJ ' 5 > 4 \u^ w M 2 ) [( S - m|) 2 + m|r|] [(m 2 + -(. + *)- m 2 w ) 2 + m 2 ^ 2 ,] 25 ' 

(56) 

?V(M*M\ = - ( g4A * \ ~ m|)(m 2 - (a + t) - m 2 ^) + m z m w T z T w 

«V hJVIS) 4 ^^2^2 ) [{s _ m | )2 + m |r|] [(m2 _ (s + t) _ ^2 + 35- 

(57) 

The total amplitude of the 3-body decays is given by the sum of all these terms, each being summed 
over all the relevant flavors and colors of the final states. The traces in equations (l43l to (I57T ) are given by, 

Tn = ^s^3m G - 3m 4 (s + 2t) + m 2 (s 2 + 8st + 6i 2 ) - s(s 2 + 2t(s + i))}, (58) 
T22 = J^( c v + c\)s\?,m & - 3m A (s + 2t) + m 2 (s 2 + 6t 2 + 8st) - s(s 2 + 2t(s + t))}, (59) 
T 33 = ^{(4 + c|)(m 2 - s)(?n 2 (2s + t) - t(s + t)) - 2c v c A m 2 s(m 2 - s - 2i)}, (60) 
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Tia = ^2 c vs{3m 6 - 3m 4 (s + 2t) + m 2 (s 2 + 8st + 6t 2 ) - s(s 2 + 2t(s + t))}, (61) 



32 r ~i 

Ti 3 = s<^ cy(3m 4 - 2m 2 (s + t) - s 2 + 2s£ + 2t 2 ) - Cj4 (3m 2 - s){m 2 -s-2t)\, (62) 

3m I. J 



32 r 

T 23 = ^ s {(4 + 4)( 3 ™ - 2m 2 (s + t) - s 2 + 2i 2 + 2st) 
- 2c v c A (3m 4 - 2m 2 (2s + 3t) + s 2 + 2si)}, 



T 44 = ^-^s|3m 6 - 3m4 ( s + 2t ) + m2 (« 2 + 8s * + 6 * 2 ) - s ( s2 + 2st + 2i2 )}> 

128 ( ~i 

T 55 = ^— ^ |m 4 (3s + t) - m 2 (s + t)(3s + t) + st(s + t)|, 

T 45 = — shm 4 - m 2 (3s + At) + t(2a + t)\, 
3m I J 

r i4 = ^2 s { 3m6 - ™ 4 ( 3s + 7t ) + "i 2 t(4s + 5t) - t 2 (s + t)}, 

T 24 = ~ CA)s{3m 6 - m 4 (3s + 7t) + m 2 t(4s + 5i) - t 2 (s + 

T 34 = ( c v - c j4 ){m 5 (3s + t) - m 3 (3s 2 + 6st + 2t 2 ) + mt(s + t)(2s + 

Ti 5 = — sf3m 4 -m 2 (3s + 4t)+t(2s + t)|, 
3m I- J 

r 25 = P~(cv - c A )s|3m 4 - m 2 (3s + At) + t(2s + t)\, 
Am L J 



'35 



64 
3m 2 



(cy - CA)|(m 2 - s)(m 2 s + (m 2 - t)(s + £)) + m 2 s(m 2 - s - 2i)j. 



(63) 

(64) 
(65) 
(66) 
(67) 
(68) 
(69) 
(70) 
(71) 
(72) 



Finally those amplitudes are applied to the golden rule for decays in order to obtain the partial and total 
decay rate for each process G — > 2 + 3 + 4 



dT = \M\' 



S 



2m 



3/2 



d 3 p 2 d 3 p 3 d 3 p4 
(2vr) 3 2£ 2 (27r) 3 2£ 3 (2tt) 3 2£ 4 



(2tt) 4 <5 4 ( Pi -p 2 - P3 - P4 ) . (73) 
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